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ABSTRACT 

We report the detection of the CO J = 1 — emission hne in three near-infrared selected star- 
forming galaxies at z ~ 1.5 with the Very Large Array (VLA) and the Green Bank telescope (GBT). 
These observations directly trace the bulk of molecular gas in these galaxies. We find H2 gas masses 
of 8.3 ± 1.9 X 10^° Mq, 5.6 ± 1.4 x 10^° Mq and 1.23 ± 0.34 x 10" Mq for BzK-4171, BzK-21000 
and BzK-16000, respectively, assuming a conversion aco = 3.6 Mq (K km s~^ pc^)~^. We combined 
our observations with previous CO 2 — 1 detections of these galaxies to study the properties of their 
molecular gas. We find brightness temperature ratios between the CO 2 — 1 and CO 1-0 emission 
lines of 0.80toJ2> l-22J:o:36 and 0.4i;^i^ for BzK-4171, BzK-21000 and BzK-16000, respectively. At 
the depth of our observations it is not possible to discern between thermodynamic equilibrium or 
sub-thermal excitation of the molecular gas at J = 2. However, the low temperature ratio found for 
BzK-16000 suggests sub-thermal excitation of CO already at J = 2. For BzK-21000, a Large Velocity 
Gradient model of its CO emission confirms previous results of the low-excitation of the molecular 
gas at J = 3. From a stacked map of the CO 1-0 images, we measure a CO 2 — 1 to CO 1-0 



brightness temperature ratio of 0.92tQ^g. This suggests that, on average, the gas in these galaxies 
is thermalized up to J = 2, has star-formation efficiencies of ^ 100 Lq (K km s~^ pc^) 
consumption timescales of ^ 0.4 Gyr, unlike SMGs and QSOs at high redshifts. 
Subject headings: galaxies: evolution — galaxies: formation — cosmology: observations - 
starburst — galaxies: high-redshift 



'^ and gas 
- galaxies: 



1. INTRODUCTION 

A remarkable step in our understanding of galaxy for- 
mation has come from the determination of the unob- 
scured history of the star formation rate (SFR) den- 
sity over a wide range of redshift. The SFR history 
of the Universe has a peak from z = 3 to 1 , and 
presents a steady decline at z < 1 (e.g. iLillv et al.l 
119961: iMadau et all 119961: iSteidel et afl 119990 . Appar- 
ently, a significant contribution to this peak is pro- 
vided by the vigorous bursts of star formation trig- 
gered by galaxy mergers and interactions. In these 
cases, the star formation is very efficient, showing high 
SFRs and high gas surface densities on scales of sev- 
eral kpc, similar to wha t is observed in distant su bmil- 
limeter galaxies (SMGs; iTacconi et al."2006', '200 8D and 
quasars (QSOs; Walter et al. 2004, 2009; Ricchcrs eFaLl 
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120081120091 ). However, an important fraction of the SFR 
density appears to come from normal star-forming disk 
galax ies (|Bell et al.ll2005l: lElbaz et al.ll2007l: IGenzel et al.l 
120081 ). Here, continuous flows of cold gas from the inter- 
galactic medium may provide the necessary fuel for star 
formation (Kcrcs ct al. 2005; Dckcl et al. 2009). 

Studies of the properties of the molecular gas 
(e.g., CO) in galaxies during the main epoch of 
galaxy formation {z — 1 — 3) have mostly focused 
on SMGs and QSOs due to the limited capabilities 
of m illimeter and radio telescopes and instruments 
(e.g., [S olomon fc Vanden Bout 2005; 'Wcifi ct al." "2005; 
Riechers et al. 2006; CariUi ct al. 2007; WciB ct al. 2007 



Aravena et alll2008l: ICoppin et all 120081: ICarilli et al] 
20101 ). These IR-luminous objects were found to have 
large reservoirs of molecular gas (M(H2) '^ lO^^^^^A/©) 
that sustain SFRs of - 500 - 1000 Mq yv~^ for < 100 
Myr. Similar to what is observed in local ultra-luminous 
infrared galaxies (ULIRGs), the molecular gas is in lo- 
cal thermodynamic eq uilibrium (LTE ) up t o high- J CO 
trans itions (J > 3; IRiechers et all 120061: IWeifi et al.l 
|2007( ). This is expected for warm, dense H2 molecular 
gas. 

Recently, large amounts of molecular gas, similar to 
those observed in bright high-redshift SMGs and QSOs, 
were fou nd in six relatively quiesc ent disk galaxies at 
z - 1.5 ()Daddi et al1 l2008( '20 10_a|W 4 similar galaxies 
at z r^ 1.2 and z ~ 2 (Tacconi et al.l 120101) and three 
disk galaxies at z = 0.5 (Daddi et al.ll2010bl ; Salmi et 
al. in preparation). A detailed analysis of the proper- 
ties of these near-IR selected galaxies revealed SFRs ~ 
50 - 200 Mq yr-\ with stellar masses of ~ lO^""" Mq. 
These galaxies apparently have high molecular gas frac- 
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TABLE 1 

Summary of the VLA observations 



Source R.A." Dec." Frequency Bandwidth 
(J2000) (J2000) (GHz) (MHz) 



Beam 



Pos. Ang. 



Gov. Fraction'^ 

(%) 



Rms" 
(^Jy/beam) 



BzK-4171 12 36 26.516 +62 08 35.35 46.760 

43.340 

BzK-21000 12 37 20.597 +62 22 34.60 45.710 

43.340 

BzK-16000 12 36 30.120 +62 14 28.00 45.653 

45.635"= 



100 


1.69" X 1.56" 


+09.8 


100 


2.02" X 1.67" 


-53.0 


100 


1.79" X 1.61" 


+18.2 


100 


2.03" X 1.67" 


-54.2 


X 3.125 


1.96" X 1.60" 


+48.9 


100 


1.90" X 1.62" 


+44.3 



96 
98 
76 



44 
62 
42 
72 
145 
181 



" VLA 1.4 GHz position from IMorrison et ahl l|2010l) . ^ Gentral frequency between the two 50 MHz IFs. '^ Fraction of the CO Hne covered 
by the VLA observations. '^ Image noise level measured around the source position. *= Combined continuum emission from two 50 MHz 
bandwidth channels at 45.585 GHz and 45.685 GHz. 



62" 08' 32" - 




12" 36™ 27.0^26.8-^ 26.6^ 26.4^ 26.2 
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62" 14' 24" - 



12"37'" 112^11.0^ 10.8^ 10.6^ 10.4^ 102^ 10.0^ 
Riglit Ascension (12000) 




12" 36™ 30.6*30.4** 30.2* 30.0* 29.8* 29.6* 
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Fig. 1. — VLA images of the G O 1-0 emission line f rom our sources. Contour levels are: -2, 2, 3 and 4a. The crosses indicate the VLA 
1.4 GHz position for our sources IMorrison et al.ll2010l) . 



tions , low star- formation effi ciencies (SFEs) (Dadd i et alj 
[20081 I2010at iTacconi et aLll201Qi and CO luminosity to 
gas mass conversion factors similar to that found in the 
Milky Way galaxy di sk rather than that of starbursting 
ULIRGs and SMGs (|Daddi et al.ll2010aD . This means 
that these objects have yet to convert a large fraction 
of their gas into stars, indicating galaxies right in the 
process of stellar build-up. The first CO 1 — mea- 
su rements in two BzK galaxi es were previously reported 
bv IDannerbauer et al.l (|2009l ). but only led to a stacked 
detection at ~ 3ct. Those observations, combined with 
the first detection of the CO 3 — 2 line in one o f these 
galaxies (BzK-21000), allowed IDannerbauer et~all ()2009f ) 
to perform an analysis of the excitation properties of the 
molecular gas, finding low excitation conditions, similar 
to what is seen in local disk galaxies. This suggests an 
important difference to what is found in more luminous 
objects at high-redshift, where the CO emission is in LTE 
up to high— J CO transitions. 

In this paper, we report individual detections of the 
CO J = 1 — emission line in a sample of three near-IR 
selected star-forming galaxies at z ^ 1.5 that have pre - 
viously been studied in detail bv IDaddi et al.l ()2010aD : 
BzK-4171, BzK-21000 and BzK-16000. In section [H we 
describe our CO 1-0 observations. In section [3l we 
present our results, give fluxes and luminosities. In sec- 
tion SI we summarize our results and discuss the average 
properties of the molecular gas in our galaxies. We as- 
sume a concordance cosmology with Hq = 71 km s~^ 
Mpc-\ rix = 0.73 and fi™ = 0.27. 



2. OBSERVATIONS 
2.1. Very Large Array 

We used the Very Large Array (VLA) in its C and D- 
array configuration and the Q-band receivers to observe 
the redshifted CO 1-0 emission line (i^rost = 115.271 
GHz) in three BzK galaxies at z ~ 1.5 in the GOODS 
North field. This configuration provides good spatial res- 
olution (typical beams of ~ 1 — 2"), with a primary beam 
of ~ 60". 

For two of our sources, BzK-4171 and BzK-21000, the 
observations were done in D-configuration between 2009 
November 02 and 2009 November 27 under very good 
weather conditions. A total of 35 hrs were spent ob- 
serving each source, over 5 tracks per source. Since the 
observations were done during a transition period to the 
Expanded VLA (EVLA), many antennae were malfunc- 
tioning and we were only able to reach less than half the 
sensitivity expected for the whole array (Table [1]). 

We used two channels of 50 MHz bandwidth each and 
two polarizations per channel. At ~ 45 GHz, 50 MHz 
correspond to ^ 330 km s~^ velocity coverage. For 
BzK-4171, the two channels were centered at 46.735 and 
46.785 GHz. For BzK-21000, the two channels were cen- 
tered at 45.685 GHz and 45.735 GHz. A fraction of the 
time of each track was spent observing these sources at 
43 GHz using two non-overlapping channels of 50 MHz 
bandwidth (i.e. 100 MHz bandwidth in total) in order to 
obtain a limit for the continuum emission. In both cases 
the phase tracking center was pointed about 10" south 
from the target positions. 
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TABLE 2 

Properties of the molecular gas. 






Source 


'^43GHz 

(MJy) 


■^CO 1-0 

(/^Jy) 


-'^CO 1-0 
(Jy kms-l) 


r; d 
(xlO"' K kms-i pc2) 


(xlQiOAfo) 


(-gas 

(xIQS yr) 


SFE9 
(L© (K km s-i pc2)-i) 


BzK-4171 

BzK-21000 

BzK-16000 


< 125 

< 144 
<360 


305 ± 70 
180 ± 50 
805 ± 220 


0.20 ±0.05 
0.13 ±0.03 
0.28 ±0.08 


2.31 ±0.53 
1.55 ±0.40 
3.42 ±0.94 


8.3 ±1.9 
5.6 ±1.4 
12.3 ±3.4 


0.80 
0.25 
0.81 


43 

142 
44 



" Upper limit (2(t) to the continuum emission at 43.34 GHz for BzK-4171 and BzK-21000 and at 45.6 GHz for BzK-16000. *> Flux density 
integrated over an area that matches the CO 2 — 1 measurements. '^ Velocity integrated intensity, I = f Sdv, corrected for emission that 
falls outside the range covered by the VLA band. CO luminosity. '^ H2 mass computed using a CO luminosity to gas mass conversion 
factor of 3.6 (K km s~^ pc^)~^ IDaddi et al.|[2010al') . J Gas depletion lifetimes and ^ star formation efficiencies derived using the SFR 



and the far-IR luminosities from Paddi et aU l)2010al '). and our CO 1 — measurements. 
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Fig. 2.— GBT spectra of BzK-4171 and BzK-21000. Horizontal 
marks indicate the position of the two 50 MHz channels observed 
by the VLA in each source. The line frequency obtained with the 
Gaussian fit (see text) was used to set the velocity scale. The 
dashed line shows single Gaussian fits to the spectra (see text), 
which were used to set the velocity scales. 

The observations of BzK-16000 were done in C- and 
D-configuration between 2009 July and 2009 December 
under mostly good weather conditions. A total of 48 
hours were used to observe this source. Spectral line ob- 
servations were performed using two IFs of 7 channels 
each and two polarizations. Both IFs have one chan- 
nel overlap, leading to a total of 13 independent spectral 
channels. These channels were combined into one single 
channel of 40.625 MHz bandwidth and centered at 45.653 
GHz, covering most of the expected CO 1 — emission 
line. We also performed continuum observations of this 



source by placing two 50 MHz channels at each side of 
the expected CO 1 — line. These channels were centered 
at 45.585 GHz and 45.685 GHz, respectively. The higher 
frequency channel overlaps the central 40.625 MHz band- 
width region by 18 MHz. 

In all cases, we used fast-switching calibration and 
observed the VLA calibrators J1302-I-5748 and 3C286 
(J1331-f305) for flux calibration. The AIPS software was 
used for data editing and calibration. Most of the data 
showed good phase stability; however, some time ranges 
were removed due to antennae with bad amplitudes. Fi- 
nally, we used the AIPS task IMAGR, which employs the 
CLEAN algorithm, and natural weighting to deconvolve 
the images down to residuals of ~ Ict in a box centered 
on our targets. Table [1] summarizes the relevant obser- 
vational parameters. 

2.2. Green Bank Telescope 

We also performed observations of two of our BzK 
galaxies with the Robert C. Byrd Green Bank telescope 
(GBT) during 2009 April and 2009 October to Novem- 
ber, under very good weather conditions. Typical com- 
puted opacities at 45 GHz are r ^ 0.1 — 0.2 and measured 
wind velocities were < 2 m s~^. At these frequencies the 
GBT beam is ^ 16". We observed in sub- reflector nod- 
ding mode, with a half-cycle time of 6 s. At the begin- 
ning of each run, we observed the following flux density 
calibrators for 10 min each: 3C286, 3C147 or 3C295. 
We estimate our measured fluxes densities to be accu- 
rate within ±10%. We used the source 1259+514 as the 
pointing and focus calibrator. Pointing and focus were 
checked every 30 min to 1 h depending on the stabil- 
ity and accuracy of the obtained corrections in each run. 
Pointing was stable within 5" in all runs. Typical sys- 
tem temperatures at 45.7 GHz and 46.7 GHz were in the 
range 85-110 K. We employed two IFs of 800 MHz band- 
width each and two polarizations each, with a spectral 
resolution of 390.625 kHz per channel, or ~2.6 km s~^ 
per channel. We placed the center of both IFs 200 MHz 
from each other and about 100 MHz from the line fre- 
quency each. In this way, the overlap region covered a 
band of -600 MHz. 

The data were reduced using the GBTIDL software, 
removing a few scans that presented bad channels or 
strongly distorted baselines. For each source, we aver- 
aged all scans and both IFs. The combination of both 
IFs produce a gain of about 5 — 10% in the obtained 
signal-to-noise ratio due to digitization noise. We fitted 
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Fig. 3.— PdBI CO 2 - 1 spectra (in arbitrary flux units) of BzK-4171, BzK-21000 and BzK-16000 (from Daddi et al. 2010a), shifted to 
the frequency of the CO 1-0 transition to illustrate the fraction of the line covered by the VLA observations. A Gaussian fit to the emission 
is also shown. The bars on top as well as the shaded area in each spectrum indicates the spectral line coverage of the VLA bands. The 
lower red horizontal bars in the BzK-16000 panel show the two 50 MHz bands used to measure the continuum emission for that galaxy. 




Fig. 4. — VLA CO 1-0 emission line maps overlaid on the HST Biz color images. Images are 4.2" X 4.2" in size. Contours and cross 
symbol are as in Fig. [T] 



polynomials of order 3 and subtracted them from the 
averaged spectra. This eliminates baseline structure on 
scales ^ 2 times larger than the expected linewidth of 
the CO lines. Finally, we downgraded their spectral res- 
olution to 61 and 63 km s"! (- 9.5 MHz) for BzK-21000 
and BzK-4171, respectively. 

3. RESULTS 

Figure [1] shows the VLA maps of CO 1 — line emission 
and Table [2] summarizes the results. Fitting two dimen- 
sional Gaussians to the CO images suggests marginally 
resolved sources, although the signal-to-noise is insuffi- 
cient to provide accurate estimates of source sizes. In 
order to derive the total CO fluxes of the sources, we 
fit Gaussians that were constrained by the position and 
size of the sources derived from the CO 2 — 1 obser- 
va tions at the Plateau de Bure Interferometer (PdBI) 
by iDaddi et al.l (|2010a| ) . Constraining the Gaussian fit- 
ting decreases the number of free parameters and permits 
a direct comparison of amplitudes, although it assumes 
that the CO 1 — and 2—1 emission are spatially coin- 
cident. 

We find that all sources are detected with significances 
in the range 3 — 4(t in the CO maps, with no evidence for 
emission in the 43 GHz continuum maps of BzK-4171 
and BzK-21000 nor in the 45 GHz continuum maps of 
BzK-16000. We derive 2cr hmits of 125, 144 and 360 ^iJy 
for the continuum flux densities of BzK-4171, BzK-21000 



and BzK-16000, respectively. 

Figured shows the GBT spectra. BzK-4171 and BzK- 
21000 are both marginally detected. Using a single 
Gaussian fit, we find peak flux densities of Sco i-o = 
0.30 ± 0.14 mJy and Sqo i-o = 0.22 ± 0.10 mJy, where 
the quoted errors account for a 10% uncertainty in the 
flux calibration, linewidths of 430 ± 190 km s^^ and 
480±220 km s'^, and CO redshifts of 1.465±0.003 and 
1.521 ±0.003 for BzK-4171 and BzK-21000, respectively. 
Our GBT measurements are consistent with those ob- 
tained with the VLA in terms of flux densities, and with 
those obtained with the PdBI for the CO 2 — 1 line in 
terms of linewidths and central frequencies. Given the 
lower S/N of the GBT observations, we refer to the VLA 
flux measurements in the remainder of this paper. 

FigureBlshows the PdBI CO 2—1 spectra of our sources 
(jDaddi et al.ll20103 ) scaled to the frequency of the CO 
1-0 transition. The horizontal bar and shaded region 
of the emission lines represent the velocity (and/or fre- 
quency) range covered by our VLA CO 1-0 observations. 
The velocity ranges covered are 641 km s~^, 656 km 
s^i and 266 km s'^ in BzK-4171, BzK-21000 and BzK- 
16000, respectively. This implies that we cover 96%, 98% 
and 76%, respectively, of the CO 1 — line with our ob- 
servations (Table [T|). After correcting for the fraction of 
the line that falls outside our bands, we derive integrated 
line intensities, Iqo = f Scodv, of 0.20±0.05, 0.13±0.03 
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Fig. 5.— CO excitation ladder or CO SED for BzK-21000 (solid 
black circles). The black solid, dashed and dotte d lines show repre - 
sentative LVG models 1, 2 and 3 from Danncrb auer et al.l II2009I ). 
respectively. The open stars and open circles show the CO 1—0 and 
CO 2 - 1 integrated flux densities for BzK-4171 and BzK-16000, 
scaled to the CO 2 — 1 emission of BzK-21000 for comparison. The 
open triangle, open square and asterisk symbols illustrate the CO 
SED of the Milky Way, the Antennae overlap region and the spiral 
galaxy NGC6946, normalized to the CO 2 — 1 emission of BzK- 
21000. The LVG models that reproduce the data for these objects 
are shown as gray solid lines. A small horizontal shift has been 
applied to the data points to enhance the visibility. 

and 0.28 ± 0.08 Jy km s^^ in BzK-4171, BzK-21000 and 
BzK-16000, respectively. We use the integrated line in- 
tensities to compute the CO line luminosities through 



L'co = 3.25 X 10^(1 



^^L^obt^CO, 



where Dl is the 



luminosity distance at rcdshift z and i^obs is the observed 
frequency (^Solomon et al. 1997 ). 

A comparison of the CO 1 — maps with optical and 
the radio positions (Fig. 2]) shows a small offset of ~ 0.3" 
from the CO peak position in the case of BzK-4171. This 
is well within the range expected given the low signifi- 
cance of the detection. For the other sources, the CO 
emission is consistent (< 0.2") with the position of the 
radio and optical source. 

3.1. Gas properties 

Our observations of the CO 1-0 line emission can be 
used to directly measure the amount of molecular gas 
and star formation efficiency in the BzK galaxies. The 
estimation of the molecular gas masses is typically done 
by using a CO luminosity to gas mass conversion factor, 
aco- Comparing observations with simulations of star- 
forming disk galaxies that reproduced the observed CO 
2 — 1 line shapes and UV morphologies o f three BzK 
galaxies (two of which are in our sample), iDaddi et alj 



(|2010aD estimated aco = 3.6±0.8 M© (K km s'^ pc^)" 



close to the Galactic disk a. Assuming this value of aco, 
we compute the H2 masses given in Table [2] 

The star formation efHciency of galaxies is defined as 
the ratio between the IR luminosity and the CO 1 — 
luminosity, SFE= L\-R/L'nr\- Using th e values for the 
IR luminosity from IDaddi et al.l ()2010aD . we find SFE= 
43, 142 and 44 L© (K km s'^ pc^)"! for BzK-4171, 
BzK-21000 and BzK-16000, respectively. We can also 
compute the time in which the gas would be consumed if 
the current SFR remains const ant, T^as = M(H 2)/SFR. 
Using the SFRs found bv .Daddi et all (|2010aD , we find 



gas consumption lifetimes of ^ 0.2 — 0.8 Gyr (see Ta- 
ble [It . These value s are comparable to those found by 
IDaddi et al.l ()2010aD and are similar those found fo r com- 
parable galaxies at redshifts ~ 1— 3 (|Tacconi et al.ir2010.) . 

3.2. Excitation of the molecular gas 

In this section we study the properties of the molec- 
ular gas in our BzK galaxies. All of our sources have 
previously been detected in the CO 2 — 1 emission line. 
However, only BzK -21000 has been detected in the CO 
3 — 2 line (Danncrb auer et"ani2009D . 

Figure [5] shows the velocity integrated CO line fiuxes of 
BzK-21000 as a function of rotational quantum number, 
J. Also shown are the CO 1 — and CO 2 — 1 integrated 
line fluxes of BzK-4171 and BzK-16000, normalized to 
the CO 2 — 1 intensity of BzK-21000 for comparison, 
and the normalize d CO inte nsities for the inner disk of 
the Milky_Way (Fix sen et al.i:1999) . the Antenna overlap 



region (IZhu ct al. 20031) and the spiral galaxy NGC6946 
(jBavet et al., i2006'). 

From the velocity integrated line fluxes, we compute 
the brightness temperature line ratios as r2i = T21/T10 = 
{l2i)/Iio) X (^^10/^21)'^, where T21 and Tio are the bright- 
ness temperatures, /21 and /lo are the integrated fluxes, 
and 1/21 and z^io are the observed frequencies of the CO 
2 — 1 and 1 — emission lines, respectively. For LTE, 
we expect this ratio to be r2i = 1. In our sources, we 
measure brightness temperature line ratios of 0.80lg'22, 

1.22l°J^ and 0.4lj:o;i3 for BzK-4171, BzK-21000 and 
BzK-16000, respectively. 

At the significance of our detections, it is not possible 
to conclude whether the emission in the individual BzK 
galaxies is thermalized or not up to J = 2. In BzK- 
4171 and BzK-21000 the values of r2i suggest thermal 
equilibrium, while in BzK-16000, the low value for r2i 
suggests that the CO emission is sub-thermal at J = 2. 
However, this needs to be confirmed with deeper obser- 
vations. For BzK-21000, the brightness temperature ra- 
tio between CO 3 - 2 and CO 1 - is rgi = 0M°:1%q, 
compatibl e with the previous result s of r^i ^ 0.5 in this 
galaxy bv iDannerbauer et al.l (|2009( ). 

For BzK-21000, we compute a Large Velocity Gra- 
dient (LVG) model using our new CO 1-0 measure- 
ment. We employ a single component LVG model that 
assu mes spherical geometry. We use the collision rates 
from lFlowed ()2001f ) with an ortho-para H2 ratio of 3 and a 
CO abundance per velocity gradient [CO]/{dv/dr)—lO~^ 
pc (km s-i)-i (e.g., |Weifiet al. 2005, 2007). Values that 
resemble the CO emission are in the range Tkin = 20—150 
K and n(H2) = 400 - 2500 cm'^. Since we do not 
have sufficient constraints to fit the data with a specific 
model, we thus discuss the three repr esentative models 
presented bv IDannerba uer et al.l (120091 ) . as shown in Fig. 
El These models have Tkin = 25 K, n(H2) = 1300 cm"^ 
and a cloud filling factor of ~2% (Model 1); Tki„ = 90 
K, n(H2) ~ 600 cm^'^ and a similar cloud filling fac- 
tor (Model 2); and Tki„ = 10 K, n(H2) = 2500 cm-"^ 
and a filling factor of 8% (Model 3). We see that all 
these models can reasonably reproduce the data for BzK- 
21000. This is expected as our CO 1 — flux measurement 
confir ms the LVG-based prediction of lDannerbauer et al.l 
(|2009f) . Jco 1-0 ~ 0.15 Jy km s'^. Therefore, we ver- 
ify that models with high Tkm and/or n(H2), as seen in 
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Fig. 6. — Combined CO 1 — emission from the three BzK galax- 
ies. Contours levels are; -2, 2, 3, 4, 5 and 6cr. The cross indicates 
the reference VLA position. 

ULIRGs or high-redshift QSOs, arc unlike ly for t he gas in 
this galaxy, as suggested by Danncrbaue r et al.l (|200&) . 

4. SUMMARY AND DISCUSSION 

We have detected the CO 1 — emission line in three 
massive star- forming galaxies at z '^ 1.5. Our obser- 
vations allow us to carry out a direct comparison with 
studies of galaxies at low-redshift, and hence put unique 
constraints on the properties of the molecular gas in disk 
galaxies at high-redshift. 

From our molecular gas excitation analysis for BzK- 
21000, we confirm the Dannerbaucr et al. (2009) results 
that the gas in this galaxy has low excitation conditions 
at J = 3 (Fig. [S])- The CO ladder for this galaxy seems 
to be thermalized up to J = 2 and it is sub-thermal 
at J = 3. The CO ladder is simil ar to that found for 
local disks, as in NGC6946 (Fig. [5] iMauersberger et aD 
Il999f ). However, at the significance of this detection, the 
CO ladder is also consistent with the Milky Way galaxy, 
for which the emission appears to be non-LTE already at 
J = 2. 

A similar case is that of BzK-4171, where the relatively 
large uncertainty in the temperature ratio, r2i, does not 
allow us to differentiate between thermal or sub-thermal 
molecular gas up to J = 2. Measurements of the CO 3 — 2 
transition are necessary to check whether or not higher 
order CO transitions are in LTE in this galaxy. Inter- 
estingly, in the case of BzK-16000, our measurements 
suggest that the CO emission is sub-thermal {r2i < 1) 
already at J = 2, resembling what appears t o be the case 
for th e inner disk of the Milky Way galaxy (|Fixsen et al.l 
|1999[) . Our estimate of rai ~ 0.4 differs from the LTE at 
the w 2(7 level. Measurements of higher— J transitions 
will help to study this particular galaxy in more detail. 

Base d on the previous CO 1 — observations of BzK- 
210 00 (iDannerbauer et all [2001. iDaddi et al.l (|2010a| ) 
and lTacconi et al.l (J2010D assumed r^i = 0.5 to convert 
their CO 3 — 2 luminosities into CO 1-0 luminosities. 
For this galaxy, we obtain r^i ^ 0.6, which validates, 
within 1(7, their assumption. 

For a solid measurement of the average CO 1 — line 



emission and the average molecular gas properties of our 
BzK galaxies, we extracted 10" x 10" cutouts centered at 
the VLA peak positions and produced a noise-weighted 
average CO 1-0 map of their combined emission (Fig. 
[S]). We obtain a detection with a flux density of 225 ± 
36 /iJy (6.25a). Stacking with respect to the CO position 
yields a flux density of 220 ± 40 /xJy (5.5cr), consistent 
within the uncertainties (< la) with the flux measured 
by stacking with r espect to the VLA po sition. From the 
values presented in IDaddi et al.l (|2010a| l for the CO 2 — 1 
emission, we compute a noise- weighted average CO 2 — 1 
to CO 1 — brightness temperature ratio of 0.92lg^g, 
suggesting that the CO SED from these galaxies is ther- 
malized up to J = 2. This ratio is consistent with what 
is found in local disk galaxies, which have lo wer CO lu- 
mino s ities than o ur BzK galaxies (Mauersb erger et al.l 
119991 : iBavet et al.l 12006. and references therein), but 
it is also similar to other galaxy populations with CO 
luminosities comparabl e to that of our BzK galaxies, 
such as local ULIRG s (iWeifi et al.' '2005"; 'Giistcn^eiaLl 
|2006!;l Hitschfeld et al.|[20 08: Bavct et al. 2006, 2009) , and 
high-redshift SMGs and QSOs (Solomon Sz Vandcn Boutl 
200^ Weifi et al. 2005 : Riechers et al. 2006: Wcifi et aU 
20071 ). Thus, our comparison based only on CO 2 — 1 to 
CO 1-0 ratios is not yet sufficient to decide whether the 
excitation conditions in our BzK galaxies are, on average, 
different or not fro m those o f other galaxy populations. 

We note that Da ddi et al.l (2010a) assumed r2i = 0.86 
to convert the CO 2 — 1 luminosities into CO 1-0 lumi- 
nosities for their sample o f 6 BzK galaxies, following the 
IDannerbauer et all (|2009[ ) results on BzK-21000. This 
assumption is supported by our measurement from the 
average stacked CO map. 

From our stacked CO 1-0 measurement, we find on 
average a SEE « 100 (K km s^^ pc^)"^ This, in- 
cluding individual values, is well below the average for 
high-redshift SMGs, 560 ± 210 L^;, ( K km s'^ pc^)'^ 
(|Greve et al.ll2005l: iTacconi et alll2006 ) and similar to the 



values found for local spiral galaxies ( Boselli et al.ll2002l : 
iLerov et al. 20(M) and for o ther disk galaxies at high- 
redshift (jTacconi et al.l[2010f ). However, it is marginally 
consistent with the average value f ound in local ULIRGs , 
^ 225 Lm (K km s'^ pc^)"! (ISolomon et all fl997l : 
lYao et al] l2003f) . We obtain an average gas depletion 
lifetime of ~ 0.4 Gyr for the galaxies in our sample, 
whic h is higher than wh at is found in SMGs, < 0.05 
Gyr (ITacconi et al.ll2006f l. 

Observations of disk galaxies at high-redshift (includ- 
ing this work) have provided detections only for the 
lower-J CO transitions, and thus they probe the lower 
end of the CO SED. Recent studies of high-redshift 
SMGs and QSOs have indicated that two gas compo- 
nents can be present: a cold, low densi ty gas component 
and a warni dens e component (e.g., iWeifi et 311120071 : 
iCarilli et al.ri201Cil ). If indeed such is the case for star- 
forming BzK galaxies, observations of higher- J CO lines 
are crucial. High-definition CO multi-transition studies 
over larger samples are necessary. Here, the EVLA and 
Atacama Large Millimeter Array (ALMA) will play a 
fundamental role, disentangling the molecular gas distri- 
bution of the cold molecular gas traced by CO 1-0 and 
the warm gas probed by the higher-J transitions. 
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